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Abstract. We present a new chemodynamical code based on the adaptive mesh 
refinement code ramses. The new code uses Eulerian hydrodynamics and N-body dy- 
namics in a cosmological framework to trace the production and advection of several 
chemical species. It is the first such code to follow the self-consistent evolution of 
chemical elements in cosmological volumes while maintaining sub-kiloparsec resolu- 
tion. The code will be used to simulate disk galaxies and explore the influence of 
chemical evolution models and star formation on galactic abundance ratios. 



There are numerous cosmological codes on the market but studies of chemical 
evolution (CE) are limited to smoothed particle hydrodynamics with a dearth of Eu- 
lerian implementations that include a detailed chemical evolution model. With the 
aim of providing a complementary approach to existing CE codes we present a fully 
cosmological, CE code with an adaptive mesh refinement hydrodynamics scheme that 
traces the formation and subsequent evolution of H, He , C, N, O, Ne, Mg, Si and 
Fe. The model presented here uses a iKroupa et al.l (1 19931) initi al mas s fracti on (IMF) 
and a SNIa delayed ti me distribution inspired by iKawata & Gibsonl (120031). Stellar 



lifetim es are taken from Kodama & Ari motol (Il997h . SNIa yields from Iwamot o et al. 



( 1999 ), SNII yields from Wooslev & Weaver! d 19951) and AGB stellar wind yields from 
van den Hoek & Groenewegenl d 19971) . 



Results 

Our code is applied to cosmological simulations of disk galaxies achieving a resolution 
of 436 pc. The galaxy presented is a field spiral galaxy of total mass 6.8xl0 11 M o . This 
galaxy will be used as a fiducial model for a series of test runs using different IMFs and 
SNIa models to ascertain the influence that each of these ingredients has on the abun- 
dance ratios, metallicity gradients, morphology and kinematics of the galaxy. A degree 
of success is achieved in fitting observations of the Milky Way disk (however we stress 
that this galaxy is by no means a Milky Way clone) but most parameter combinations 
are too a-rich at the high metallicity end of the distribution. It is believed that this can 
be traced to the relatively low SNIa rate and future runs will explore this in more detail. 



Summary 

We present the first of what will become a suite of Eulerian cosmological disk galaxy 
simulations with CE and sub-kpc resolution. The first simulations show a good agree- 
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Figure 1. Star formation and SN rates (left panel). The z=0 SN rates are 
SNII=0.763SJVhM and SNIa=0.096SjVt <M (SN per century per 1O 1O M ) and com- 
pare favourably with those observed by iMannucci et al] d2008l) of 0.65*!!'j|S NuM 
(SNII) and O.U^S NuM for (SNIa) in field Sbc/d galaxies. An example of 
[O/Fe] vs. [Fe/H] (right panel) for t he particular CE mo del used in this run. Ob- 
serva tional data is plotted in grey, ICavrel et alj d2004l) (crosses), F/G/K dwarfs 
from iRamfrez et alj d2007l) (diamonds) and solar neighborhood F/G stars from 
Bens by et alJ (l2004h (triangles). 



ment of the SN rates with observations and an improved rotation curve with respect to 
their counterparts created with the standard feedback mechanism. A great deal of vari- 
ation in abundance ratios is seen under changes in initial mass function slope, upper 
mass limit and SNIa delayed time distribution. Future work will explore a full range in 
parameter space to constrain the CE of disk galaxies. 
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